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The aim of the paper is to present numerical results supporting the presence of conformal invari¬ 
ance in three dimensional statistical mechanics models at criticality and to elucidate the geometric 
aspects of universality. As a case study we study three dimensional percolation at criticality in 
bounded domains. Both on discrete and continuous models of critical percolation, we test by nu¬ 
merical experiments the invariance of quantities in finite domains under conformal transformations 
focusing on crossing probabilities. Our results show clear evidence of the onset of conformal invari¬ 
ance in finite realizations especially for the continuum percolation models. Finally we propose a 
simple analytical function approximating the crossing probability among two spherical caps on the 
surface of a sphere and confront it with the numerical results. 


I. INTRODUCTION 

A major fact in the theory of critical phenomena is that the scale invariance exhibited by systems at 
criticality [T] when combined with a local action may give rise to invariance under the larger group of 
conformal transformation which locally acts as scale transformation [5]. Although being scale invariant 
and having a local action is not in general equivalent to being conformal invariant, since there are scale 
invariant local models which are not conformal invariant and non-local models possessing conformal 
invariance mzi, it is widely believed to hold true for many models of physical interest. 

The conformal group in d spatial dimensions (for d ^ 2) has a number of independent generators equal 
to ^{d-\-l){d-\-2). For the special case d = 2 the conformal group turns out to be infinitely dimensional 
This infinite dimensionality is indeed at the root of the success of conformal field theories in the study of 
two dimensional critical models. For higher dimensions the exploitation of conformal symmetry (if present) 
leads nonetheless to important simplifications in the study of critical models. Even though for a long lapse 
of time traditional techniques not explicitly incorporating conformal invariance were used for the study 
of critical phenomena (such as the e-expansion [8]), the last years witnessed the development in three 
dimensions of the so called conformal bootstrap program which, fully exploiting conformal invariance, has 
lead to high-precision nonperturbative predictions for anomalous scaling dimensions of the Ising model 
in three dimensions [S] and in fractional dimensions m- The role of conformal symmetry in the theory 
of critical phenomena can thus be hardly underestimated. 

Dating back to the Polyakov’s hypothesis of the conformal invariance of critical fluctuations m and 
to the derivation under broad conditions in d = 2 of (conformal) invariance under local changes of 
length scale from the (scale) invariance under rigid length scale changes [12], a significant amount of 
work was by then devoted to show with mostly analytical different techniques the possible existence 
of conformal symmetry at criticality in higher dimensions and its relation with the scale invariance. It 
was soon pointed out that unitarity is a key ingredient, with supersymmetric models in d = 4 being an 
important example in which the conformal invariance was explicitly shown m- Always in d = 4 a proof 
of conformal invariance within perturbation theory has been recently discussed nms]. In d = 3 we 
mention the recent discussion of the occurrence of conformal invariance for the Ising model at criticality 
using functional renormalization group techniques Hg. 

For two dimensional systems, complementary to analytical tools and techniques, an important role 
in elucidating the geometric aspects of universality and in testing the emergence and consequences of 
conformal invariance has been played by numerical simulations: in d = 2 conformal symmetry has been 
extensively tested numerically (see e.g. the early numerical experiments in |17| and CHI) by confronting 
numerics with theoretical predictions arising from conformal held theory. 

However, in other space dimensionlities these tests and numerical experiments are still partial. An 
active research line dates back to the paper |19| . in which a relation among correlation functions defined 
on a flat d-dimensional geometry and a suitably chosen curved geometry (the hypercylinder S‘^~^ x K) 
has been put forward, which in d = 2 amounts to the conformal mapping between the (both flat) 
plane and cylinder geometries. This relation was originally verified for the spherical model for dimension 
2 < d < 4 jig. Recent Monte Garlo numerical studies established the validity of this prediction for the 
three dimensional Ising model |2()l |2I] : in these works a continuous limit of the lattice model was used to 
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perform the simulations of the curved space and the assumption of conformal invariance, together with 
the determination of the magnetic and energy correlation lengths, was shown to be enough to correctly 
reproduce the known critical exponents. Similar numerical computations were carried out for the three 
dimensional bond percolation [22] and 0{N) models |28| . 

The aim of this paper is to provide numerical evidence of the invariance of three dimensional critical 
systems under conformal transformations of the euclidean flat space (more precisely, of the flat space 
into itself). For this reason we studied observables of three dimensional percolation at criticality, study¬ 
ing models both on lattice (site percolation) and on the continuum (overlapping spheres). We decided 
to choose percolation for several reasons: percolation is possibly the simplest model exhibiting a phase 
transition, see e.g. |23| and the recent review [2^. This simplicity has on one side allowed to mathe¬ 
matically prove many conjectured properties of percolation such as the existence of a transition in the 
thermodynamic limit |2()] or, more recently, to prove rigorously many results from CFT with help 
of of Stochastic Loewner Evolution (SLE) techniques |2Sj. On the other hand, with the aid of powerful 
algorithms |22l 150], the numerical exploration of percolation problems is especially convenient. Finally, 
it is clear that percolation is an ideal system to test conformal invariance in bounded domains, as con¬ 
vincingly showed for d = 2 |18| . since in a natural way one computes quantities defined and crucially 
dependent on the form of the boundaries, as the probability of connecting two disconnected regions of 
the boundary of the system. 

With these motivations, we numerically study the presence of conformal invariance for three dimensional 
percolation at criticality in bounded domains. We introduce and study both discrete and continuous 
models of critical percolation, we test by numerical experiments the invariance of quantities in finite 
domains under conformal transformations by computing crossing probabilities. We also put forward an 
approximate expression for the crossing among spherical caps on a sphere and check it against numerical 
experiments finding very good agreement. Our results show clear evidence of the onset of conformal 
invariance in finite realizations especially for the continuum percolation models. 


II. PERCOLATION MODELS AND CONFORMAL INVARIANCE 


In its simplest incarnation percolation consists of filling a randomly and independently with probability 
p the sites of a graph and checking for the existence of large clusters. As a testbed for the verification of 
conformal invariance we will consider a classical problem in percolation theory: the existence of a cluster 
connecting two disjoint regions uji and uj 2 of the boundary of a bounded domain fl. The probability of the 
occurrence of this event is known (and under some assumptions mathematically proven |26p to tend, in 
the thermodynamic limit, to a nontrivial (i.e. different from 0 or 1) value ttx if the occupation probability 
p is set to the critical probability pc- Before introducing the finite realizations considered in this work (this 
will be done in Section III \ we explore in this Section the geometric consequences of conformal invariance 


by working in the thermodynamic limit. Due to universality these considerations do not depend on the 
specific model of percolation under scrutiny. In order to verify this model independence we will take two 
prototypical models of percolation: a discrete and a continuum one, i.e. site percolation and continuum 
percolation of penetrable spheres, with the notations introduced in this Section being applicable to both 
models. 

The crossing probability has been the object of classical foundational papers imiiH] for the applica¬ 
tion of conformal field theory in two dimensional critical systems. In the two dimensional case indeed 
full exploitation of the (infinite dimensional) conformal group has allowed to obtain the exact crossing 
probability for simply connected domains in the paper by Cardy HZ]. In that setting the presence of 
boundaries gets translated, in the boundary CFT language, into the insertion of boundary creating op¬ 
erators located at the points at the extremes of the one dimensional boundary. Although many similar 
formulas in different and more complicated settings have been derived m, some details of the operator 
content in percolation in two dimensions, which turns out to be a logarithmic CFT, are still matter of 
investigation |82| . 

In our three dimensional setting the domain D will be chosen as a simply connected closed domain and 
the two surface domains as intersection of the boundary dil and two closed domains and D 2 . Although 
the two domains and 172 are not strictly needed to define of the regions wi uj 2 they will come in handy 
to define unambiguously the finite lattice realization we will study numerically. 

At criticality the crossing probability should take a value ttx (17, f7i, ^ 2 ) dependent on the three domains 
17, f7i and 172. Please note that the actual value of ttx (17, f7i, 172) does not depend on the percolation 
model we are studying. This expectation has been numerically validated in widely different models in the 
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two dimensional models and in some selected three dimensional cases |33l IM] , 

If we perform a conformal mapping £ the three domains get mapped onto the three domains £(0), 
£(ni) and £(r 22 )- The conformal invariance hypothesis implies that: 

7T^in,n,,n2) = 7T^iG:{n),G:{n,),€in2)). (i) 

Note that in general conformal invariance would allow for dimensionful prefactors in front of Q depend¬ 
ing on the geometry and on the scaling dimension of the observable under consideration. The further 
requirement on the crossing probability to be, as expected, a finite number in the thermodynamic limit 
forces us to set them to one. In the two dimensional case, using the CFT language, the above statement 
can be expressed by saying that the boundary operators we insert to obtain the crossing probability are 
expected to have zero conformal weight. In three dimensions of course the interpretation of the boundary 
condition as insertion of boundary operators is lacking but the expectation of the crossing probability to 
have vanishing scaling dimension still holds. 

Another observation is in order the above statements on implication of conformal invariance crucially 
depend on the choice of the boundary conditions chosen. Indeed different boundary conditions can flow 
to |351136] renormalization group fixed points differing from the one encountered in the bulk systems. A 
study based on e-expansion of the effect of boundary conditions on the critical properties of percolation 
models in semi-infinite geometries has been carried out in m- The boundary conditions we have chosen 
i.e. occupied for points in uji and free for the remaining part of the boundary are expected, in analogy to 
what observed in two dimension, to the so called ordinary fixed point. 

We recall that in three dimensional space the conformal group is 10-dimensional [2] and is generated 
by the following generators (acting on a three dimensional vector x): 

• 3 translations: 9^ 

• 1 dilatation: 

• 3 rotations: x^^d^, — x^d^ 

• 3 special conformal transformations: 'Y^^j(2x^Xvdi, — x^x^d^) 

where x^ denote the components of x. Indeed the special conformal transformations are most interesting 
ones and we will concentrate on them. 

We will consider different starting geometries and then examine the problem obtained by applying 
compositions of the following conformal transformations: 


, ^ (1 - v'^){x^ + v^) + u^(x -k v) ■ (x -f v) 

^ 1 -I- 2x • V -I- x^v'^ 

where the vector v rules the magnitude and direction of the transformation. We will concentrate on the two 
transformations Ca, and having vectors v = e{l, 0 , 0 } and v = e{ 0 , 0 , 1 } respectively and compositions 
of them. The parameter e ruling the magnitude of the transformation (when e = 0 transformation 
Q reduces to the identity) will be set conventionally to 0.2. These transformations, which are finite 
exponentiations of the generators —d^ -I- '^^{2x^x^di, — Xi,x^d^), i.e. combination of translations and 
special conformal transforms, have the property of mapping the unit sphere onto itself. In table|T]we write 
the transformation which will be considered together with the roman numeral used to denote it. 

We will examine the following percolation geometries: 

i) If we take the the domain to be the unit sphere x^ + < 1 and the domains defining the 

boundary domains f = 1 , 2 to be spheres intersecting with at right angles; thus the problem reduces 
to the calculation of crossing among two (disjoint) spherical caps on a sphere. The starting geometry 
(I in our numbering scheme IB is defined by two caps enclosed by parallels with azimuthal angles 9 and 
TT — 0 as depicted in Figure [bin panel a together with a geometry obtained by applying the conformal 
transformation VIII in tableU 

By use of conformal transformations it is easy to see that the crossing probability will depend only on 
one anharmonic ratio a. a is easily evaluated by mapping the surface of the sphere by a stereographic 
projection onto the (extended) plane and having the caps uji and uj 2 mapped onto two circles with center 
in the origin (zi) and the point at infinity (Z 4 ). Scale invariance allows us now to set the radius of wj to 



4 


numbering 

transformation 

I 

1 (identity) 

II 


III 


IV 

CzC. 

V 

C.Cz 

VI 

CzCz 

VII 


VIII 


IX 



TABLE I: Definition of the considered conformai transformations. 



FIG. 1: Geometries related by conformal transformations. The domain Q is the sphere centered in the origin while 
the domains Di are the the ones intersecting D which have, when necessary, been trimmed to ht the image. The 
domains w; are the the darker portions of D. Panel a refers to the cap geometry (i) with 9 — 0.457r/2 and when it 
subject to transformation VIII in table [I] it is mapped onto the geometry shown in panel b. Please note that the 
sphere Q.2 gets mapped onto the exterior of a sphere under the action of the conformal transformation. Panel c 
is the octant geometry (ii) and when acted on with transformation V in table |I] it is mapped onto the geometry 
shown in panel d. 
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be one ( 22 ). The only free parameter left is the radius of which we call Z3. We define the a to be the 
anharmonic ratio 


(Zs - Zi){z4 - Z2) \2J 


(3) 


written with its relation to the angle 0. In this geometric setting conformal invariance implies that ttx 
depends only on a. 

a) The second geometry we will consider is the crossing probability among two spherical triangles 
on the unit sphere. The domains Oi and U ,2 will be chosen as the octants x,y,z > 0 and x,y,z < 0 
respectively. The starting geometry is depicted in panel c of figure In panel d of figure we show the 
octant geometry when transformed with transformation V of Table^ 


III. PERCOLATION MODELS 

In order to test the invariance of crossing probabilities we will consider two types of percolation models, 
a discrete and a continuum one. 


A. Discrete model 

Take a (simple) cubic lattice L of lattice constant A i.e. L = \{1? + {1/2,1/2,1/2}). The shift by the 
vector {1/2,1/2,1/2} has been introduced in order to grant better scaling properties of our finite lattice 
realization. We have specialized to the site percolation problem. The sites will be filled with a probability 
p. Such a model has been one of the most widely studied among the percolation models in three spatial 
dimension. As known the percolation transition for p = pc where Pc = 0.31160768(15) |3S]. As for the 
realization of the geometries introduced in Section [T^ are concerned in the present model a point will 
be considered inside the domain if it belongs to L n 17 while the points belonging to oji are the ones in 
L n (17i \ 17) having at least a nearest neighbor in the domain. In panel a of figure]^ a 2d representation 
of the above model. 


B. Continuum model 

The continuum percolation model we have considered is the percolation of penetrable spheres [55] 
namely we extract the centers of the spheres of radius p uniformly from the interior of 17 (Poisson point 
process). A sphere will belong to the uji if the distance of its center c to 17 is less or equal than p and c 
is an element of 17^. Please see panel b of figure 

The existence of a percolating cluster depends on the so called filling factor j] defined as the mean 
number of objects n times the ratio of the volume of the filling spheres to the total volume which in our 
case reads rj = np^. The value pc for which the transition occurs is known from numerical experiments 
to be Pc = 0.34189(2) [30] and pc = 0.341888(3) jH] which is to date the most precise estimate of the 
critical filling. 


IV. NUMERICAL SIMULATIONS 

The discrete percolation experiments have been carried out on a simple cubic lattice A of decreasing 
lattice spacing, namely 1/A = 8,16, 32, 64,128 in order to perform a scaling analysis. We have employed 
the Newman-Ziff Hi |3i algorithm to calculate the crossing probability. We have thus obtained the 
crossing probability ttx (n) as the lattice is filled (incrementally) with n sites. From the above probabilities 
at fixed filling we can reconstruct the crossing probability ttx (p) in terms of any occupation probability 
p via the following binomial convolution: 


N 

(p) = (n)p”(l - p)^~ 


n—0 



(4) 
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FIG. 2: 2d cartoons of the percolation models we have considered in this work. Panel a refers to the discrete 
percolation model whilst panel b to the continuum one. The shaded circles in a represent occupied sites, while 
in b we depict the randomly occupied circles (spheres in our 3d model). The sites (in a) and circles (in b) with 
horizontal and vertical shadings (red and blue in color) represent cluster connected to the boundary domains cJi 
and 012 respectively. In both cases no cluster connecting cji and a;2 has been established. 


For the continuum percolation experiments we used spheres of radius p. We have considered the sizes of p 
to bel !p = 2^, 2^-^, 2^,..., 2®-^, 2^ with centers extracted uniformly from the volume under consideration. 

In this case we resorted to the continuum version of the Newman-Ziff algorithm detailed in [21]. In this 
case the measurement of crossing probabilities with a given number n of spheres ttx (n) allows to obtain 
the crossing probability for any filling factor ry by a poissonian convolution: 


^x(p) 


= 


^7rx(n) 

n=0 


(n/p^ 

n! 


(5) 


The number of realizations used for our simulation consists of up to 10^ samples allowing us to obtain 
the desired statistical uncertainties. 

For the random number generation we resorted to the, well tested for percolation simulations |42l I43| . 
Mersenne twister MT19937 by Matsumoto and Nishimura |44| . 


V. RESULTS 
A. Site percolation 

In figure (left panel) we depict the crossing probability as a function of the inverse lattice spacing for 
the cap geometry for some of the values of 0 considered and for the octant geometry. These have been 
obtained by setting the occupation probability p to the best known value Pc = 0.31160768(15) |38| . 

As we can see groups of curves corresponding to the geometries related by a conformal transformation 
converge to the same value as the lattice realization becomes finer. Among the various mapped geometries 
we observe that VIII and IX are the ones differing the most from the others. This is to be expected since 
they are subject to more severe deformation such that one of the caps becomes very small (cfr. panel b 
of Figure [preferring to a geometry obtained by applying transformation VIII). When this happens even 
the finest simulated finite lattice realizations have stronger finite size effects. 

In order to quantitatively assess the convergence to the same value in [p (right panel) we plot the 
variance of ttx for the nine geometries considered for different values of 0 which decay to zero, signalling 
the presence of the conformal symmetry in the thermodynamic limit. 
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FIG. 3: (Left panel) Crossing probabilities at as a function of the inverse lattice spacing 1/A. The different 
colors refer to geometries obtained by applying the transformations in table [I] The distinct bunches of curves refer 
to different values of 0 and for the octant geometry. Statistical errors are smaller than the symbols. (Right panel) 
Variance of the crossing probabilities as a function of the system size for different values of 6. The lines are just 
guide to the eyes. 


By extrapolating these finite size results to the thermodynamic limit allows us to obtain a numerical 
estimate of the crossing probabilities for any two spherical cap on the sphere. The crossing probability in 
terms of the anharmonic ratio a is depicted in In this figure we also plot the function 

TTx (a) = tanh ^tan —^ (6) 

which appears to describe very well the values extrapolated in the thermodynamic limit. Indeed fitting 
the extrapolated values the function tanh (a tan with a free parameter gives a value of a = 0.989(6) 
very close to 1. 

The function ttx we have numerically estimated and for which we have given an analytic approximation 
TTx (a) is somewhat analogous to the Cardy’s formula m although in the two dimensional case, due to 
Riemann mapping theorem, all of the crossing probabilities among two distinct segment on the bound¬ 
ary of a connected domain can be obtained by this formula. In the three dimensional case instead the 
knowledge of the function plotted in [fallows us to calculate the crossing probability among two arbitrary 
spherical caps on the surface of a sphere. 


B. Continuum percolation 

In contrast to the discrete model the continuum one has the added benefit of having a continuously 
tunable parameter p. Moreover as we lower p the various crossing probabilities approach the thermo¬ 
dynamic limit in a much smoother way. In we plot the crossing probability as function of the filling 
fraction for decreasing filling sphere radii for 9 — 0.757r/2 for the configuration I. These curves clearly 
develop a steep profile as we approach the thermodynamic limit as shown in panel a of Figure We can 
use them to obtain reliable estimates of the critical filling fraction. In fact, by defining a size dependent 
77 c(p) (the so called cell-to-cell estimator jJS]) as the value for which the curves meet: 

T^xiVcip)) =Trx'^{Vc{p)) ( 7 ) 

we can easily obtain an estimate for pc- This analysis, shown in panel b of Figure]^ with the 9 = 0.657r/2 
cap geometry I leads us to estimate the critical p as pc = 0.341935(8). We have chosen the value 9 = 
0.657r/2 and the configuration I since, among the many geometries we have simulated, it represents “best” 
the thermodynamic limit for it has the domains wi and u }2 of the same size and the ratio of sum the total 
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FIG. 4: Crossing function for two spherical caps on a sphere in term of the anharmonic ratio a (defined in[^. The 
different curves refer to different system sizes, and to the extrapolation to the infinite size system. The dotted 
lines are just guide to the eyes. The continuous line refers to the function Tt{a) dehned in Equation 0 - 


area of the patches uii and the remaining surface is closest to 1. This estimate has to be compare with the 
ones available in literature m and im which provide values of rjc = 0.34189(2) and rjc = 0.341888(3) 
respectively. Those results were obtained with growth algorithms on large systems. In our simulation 
the lowest value of 6 {9 = 0.457r/2) simulated cap geometry resembles somehow the study of percolating 
clusters connecting two sides of a cube since the ratio of patches area to open surface area is closest to 
1/2. Performing the analysis on this configuration we obtain a value of rjc = 0.34190(1) which is consistent 
with the above cited values, showing how the determination of ? 7 c sensibly depends on the geometry for 
the sizes we have examined. We will rely on our “best” determined estimate rjc = 0.341935(8) for the 
subsequent analysis but we anticipate that the final results will not be crucially affected by the choice of 
Vc- 

The convergence to the same value for ttx for geometries related by a conformal transformation can 
seen in the left panel of Figure which is the continuum analogue of the left panel in Figure The 
curves for ttx nicely converge to the thermodynamic limit, since the effects from discretization are not 
present; on the other hand finite size effects are stronger, for the range of p simulated, since their distance 
is in general bigger than the curves for the discrete systems. Once again we observe that configurations 
generated with VIII and IX transformations are more distant from the other curves confirming that, for 
the same value of p, they are more far from the thermodynamic limit. 

The good behavior of the finite realizations allows to have an independent extrapolation for each 
configuration. This is shown in the right panel of Figure]^ where the extrapolated values of the crossing 
probabilities are shown for the set of conformal transformations examined together with the mean over the 
various configurations. This confirms the onset of conformal invariance even for a continuum percolation 
model 

For comparison we also report, in the left panel of figure the same analysis performed with the 
critical filling reported in [41] for the cap geometry with 9 — 0.457r/2. As we can see in both cases the 
different transformed geometries give the same value within 1.5cr so they are consistent with each other. 
This proves that within our error, both values of pc give the same value for the crossing probability in 
domains related by a conformal mapping although the value of ttx depends on the chosen pc- Similar 
findings are obtained for the other values of 9 examined and for the octant geometry. 

The function 7r(a) 0, introduced in the previous section, ruling the crossing probability among two 
caps in on the sphere, is depicted in the left panel of|^ As we can see the values are consistent, within 
errors, with the function obtained from the site percolation model. This provides an evidence of univer¬ 
sality of TTx in percolation. Moreover if we compare it with the analytic function ttx [^proposed in the 
previous section we find an excellent agreement with the data obtained setting pc to the value provided in 
|41Jw hile if we choose our best estimate of pc we have a less faithful representation of the numerical data 
by0 If we try to fit the above data with the more general function tanh (a tan indeed we obtain the 
estimates of a = 1.022(6) for p^ = 0.341935 (our best value) and a = 1.005(5) for p^ = 0.341888 (given in 
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FIG. 5: (Left panel)Crossing probability as a function of the filling factor for the cap geometry with 9 = 0.657r/2. 
The different lines refer to the filling spheres of increasing radius 1/p = 2^, 2^'®,..., 2®'® in the direction of the 
arrow. The errors are not showed being smaller than 2.5 • 10“'* for all the points i.e. smaller than the lines size. 
(Right panel) Finite size estimations of rjdp) as a function of 1/p. The continuous line is the best fitting function 
0.341935 — 0.687988p* ®®^®®. The two lowest values of 1/p = 2^, 2^ ® have been excluded from the fit. 
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FIG. 6 : (Left panel) Crossing probabilities ttx (L) of conformally equivalent cap geometries and octant geometry 
as the radius p of the filling spheres is reduced, cfr. left panel of Figure]^ The crossing probabilies are calculated 
at 77 = Pc- (Right panel) Extrapolated values of ttx in the thermodynamic limit for the caps and octant geometry 
with for the different conformal transformations considered H] 


my 


VI. CONCLUSIONS 

By means of numerical experiments we have assessed the invariance under conformal transformation of 
selected both discrete and continuous percolation problems in three dimensions at criticality in bounded 
domains. We also proposed an analytical function approximating with very good accuracy the crossing 
probability among two arbitrary spherical caps on a sphere. 

We hope that our work can stimulate from one side a study of percolation by conformal bootstrap 
techniques and from the other further investigation of symmetries of general percolation models. Future 
work in our opinion interesting will entail the analysis of bulk observables and the study of different 
statistical mechanics models such as 0{N) and Potts models. 
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II III IV V VI VII VIII IX 

configuration 



FIG. 7: (Left panel) Extrapolated values of ttx in the thermodynamic limit for the cap geometry with 6 = 0.457r/2 
for the different conformal transformations considered [I] The two curves refer to the values obtained setting rj 
to our best estimated value (upper curve) and the value, given in |41| . (Right panel) Comparison between the 
values extrapolated at infinity for the site percolation model and the continuum percolation model with the two 
estimates for rjc as a function of the anharmonic ratio a §. We also report the function 7r(a) defined in §. 


Note added: During final stage of this work the authors noticed on arXiv a very recent and very 
interesting paper on the numerical investigation of Ising model on spherical domains with the help of the 
insight gained from the conformal bootstrap analysis of the model |33] for the finite size scaling analysis 
of correlators. The work [JB] concentrates on observables lying in the bulk of a three dimensional Ising 
model. It would be interesting to extend the analysis of |1B] to observables living on the boundary. For 
the percolation the corresponding analysis of operators living either in the bulk or in the boundary of 
the pertinent three dimensional boundary CFT has, to the best of our knowledge, not been derived. 

Acknowledgements: We acknowledge useful discussions with Defenu N and Delfino G. The authors 
benefitted from computational resources from the Iscra C project COSY3D at CINECA Bologna Italy. 
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